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1, (Canceled) 

5 2 . (Amended) The system estimation method according to claim 

j 7 or 9 or 11 , wherein the processing section calculates the existence 
condition in accordance with a following expression: 



= ^t|<-i + ~ — ^—HjHi > o, i = o,...,fc 



(17) 



10 3 . (Amended) The system estimation method according to claim 

7 or 9 or 11 , wherein the processing section calculates the existence 
condition in accordance with a following expression: 
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-Q~i + py}>0, z = 0,--,Jfc 



here. 



o = 1 - v 2 - P H i K *,i n , / s 



(18) 



(19) 



where the forgetting factor p and the upper limit value Yf have a 
following relation: 

0 < p = 1 - x(Vf) £ 1, where y(y f ) denotes a monotonically damping 
function of y f to satisfy x(l) = 1 and x(°°) = 0- 
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4. (Canceled) 

(Amended) The system estimation method according to claim 7 or 0 
or 11 , wherein the forgetting ooofficicnt p and the upper limit 
value y# have a following relation : 

0 < p - 1 — %(V*) * 1, whore %(y±) denotes a monotone damping function 
of v* to satisfy %(1) ■ 1 and %U^-«^-Q^ 
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5. (Canceled) 
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6, (Canceled) 



7. (Amended) A system estimation method for making state 
estimation robust and optimizing a forgetting factor p 
simultaneously in an estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+i - F k x k + G k w k 
y k - H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 
which indicates a ratio of a filter error to a disturbance including 
the system noise w k and the observation noise v k and is weighted 
with the forgetting factor p is suppressed to be smaller than a 
term corresponding to a previously given upper limit value y £ , and 

the system estimation method comprises: 

a step at which a processing section inputs the upper limit 
value y f , the observation signal y k as an input of a filter and a 
value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value Yf; 

a step of executing a hyper filter at which the processing 



section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s , k by using the forgetting factor p and a gain matrix 
K k and by following expressions (20) to (22), or, expression (20) 
5 and expressions which are deleted Ji" 1 and Ji in the expressions 
(21) and (22),: 



jr #Jk «jrj k (: l l)/H M (l i l) , K^pHp-lKtR^J^J^ 
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(20) 
(21) 

(22) 



Where, 



Re* = R* + Ck£k\k-iC%, Ck — 



Ji 



[o -1 



* 1 ~ T 

> Ek\k-i = ^jfc|ife-i^|jb-i 



j, R^^B^J^, xo,o = x 0 (23) 



Q(k) denotes a J-unitary matrix, that is, satisfies 
0(k) J0H(k) T - J ff J - (Ji ® X), I denotes a unit matrix, K k ( : ,1) 
denotes a column vector of a first column of the matrix K k , 

here , 

10 x~ k | k : the estimated value of the state x k at the time k using the 

observation signals y 0 to y k , 

y k : the observation signal, 

F k : the dynamics of the system, 

K s , k : the filter gain, 
15 H k : the observation matrix, 

2~k|k: corresponding to a covariance matrix of an error of x~ k | k , 

0(k): the J-unitary matrix, and 

R e , k : an auxiliary variable. 

a step at which the processing section stores an estimated 
20 value of the state x k by the hyper filter into the storage section; 

a step at which the processing section calculates an existence 
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condition based on the upper limit value Yf and the forgetting factor 
p by the obtained observation matrix Hi or the observation matrix 
Hi and the filter gain K s#i , and 

a step at which the processing section sets the upper limit 
value to be small within a range where the existence condition is 
satisfied at each time and stores the value into the storage section, 
by decreasing the upper limit value Yf and repeating the step of 
executing the hyper Hoo filter. 

8 . (Amended) The system estimation method according to claim 
7, wherein the step of executing the hyper Hoo filter includes: 

a step at which the processing section calculates and 
2~k+i|k I/2 by using the expression (22); 

a step at which the processing section calculates the filter 
gain K s , k based on an initial condition of 2~ k | k -± and an initial 
condition of C k , and tho matrix gain by using the expressions 
(21) and (22) ; 

a step at which the processing section updates a filter 
equation of the Hoo filter of the expression (20); and 

a step at which the processing section repeatedly executes 
the step of calculating by using the expression (2230), the step 
of calculating by using the expressions (21) and (22) , and, the 
step of updating while advancing the time k. 

9. (Amended) A system estimation method for making state 
estimation robust and optimizing a forgetting factor p 
simultaneously in an estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+ i = F k x k + GkWjc 
y k a H k x k + v k 
z k = H k x k 
here. 



x k * a state vector or simply a state, 
w k : a system noise, 
v k : an observation noise, 
y k : an observation signal, 
5 z k : an output signal, 

F k : dynamics of a system, and 
G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 
which indicates a ratio of a filter error to a disturbance including 
10 the system noise w k and the observation noise v k and is weighted 
with the forgetting factor p is suppressed to be smaller than a 
term corresponding to a previously given upper limit value Yf * and 

the system estimation method comprises: 

a step at which a processing section inputs the upper limit 
15 value Yf. the observation signal y k as an input of a filter and a 
value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
20 with the upper limit value Yf; 

a step of executing a hyper ft* filter at which the processing 
section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s , k by using the forgetting factor p and a gain matrix 
25 K k and by following expressions: 
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R 



0 



**|t = £*-i|fc-i + K,j.{y k - JTfc*fc_i|jfe_i) 

1 _ 1 
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(61) 
(62) 

e(Jfc) (63) 



here, 0(k) denotes an arbitrary J -unitary matrix, and C k = 
CwiW is established, where 



p* 0 

o p4 7/ 



r i o l ~ ^? 



(23) 



here # 

x~ k | k : the estimated value of the state x k at the time k using the 

observation signals y 0 to y k , 
5 y k : the observation signal, 

K s , k : the filter gain, 

H k : the observation matrix, 

0(k): the J-unitary matrix, and 

R e , k : an auxiliary variable. 
10 a step at which the processing section stores an estimated 

value of the state x k by the hyper ft* filter into the storage section; 
a step at which the processing section calculates an existence 

condition based on the upper limit value Yf and the forgetting factor 

p by the obtained observation matrix Hi or the observation matrix 
15 Hi and the filter gain K Sfl , and 

a step at which the processing section sets the upper limit 

value to be small within a range where the existence condition is 

satisfied at each time and stores the value into the storage section, 

by decreasing the upper limit value y f and repeating the step of 
20 executing the hyper H«> filter. 



10. (Amended) The system estimation method according to 
claim 9 , wherein the step of executing the hyper Ha, filter includes : 



6 



a step at which the processing section calculates K~ k based 
on an initial condition of R e ,k+i, Rr,k+i and L\+i , and the matrix gain 
K~ k by using the expression (63); 

a step at which the processing section calculates the filter 
gain K s , k based on the initial condition and by using the expression 
(62); 

a step at which the processing section updates a filter 
equation of the filter of the expression (61); and 

a step at which the processing section repeatedly executes 
the step of calculating by using the expression s (62) and (63), 
the step of calculating by using the expression (6244:) , and, the 
step of updating while advancing the time k. 

11. - (Amended) A system estimation method for making state 
estimation robust and optimizing a forgetting factor p 
simultaneously in an estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+1 = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 
which indicates a ratio of a filter error to a disturbance including 
the system noise w k and the observation noise v k and is weighted 
with the forgetting factor p is suppressed to be smaller than a 



term corresponding to a previously given upper limit value Yf / and 
the system estimation method comprises: 

a step at which a processing section inputs the upper limit 
value Yf/ the observation signal y k as an input of a filter and a 
value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value Yf; 

a step of executing a hyper H«> filter at which the processing 
section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s ,k by using the forgetting factor p and a gain matrix 
K"k and by following expressions: 



lf i j fc = p*ir*(: l l)/i^(l t l) 



Kk+l 




0 


0 







Ljt+i = p *L k 



-p-iL k R-lL k C M 

1 W " 



0 

K k 



i*r,*+l = Rrt " L k C M R €%k Ck+\L k 

Where , 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



, H k+l = [u k+lU (k + l-N)] = [u{k + l)u k ], jf, = [«(!), 0.....0] 



P o 



0 -"~ 2 



L 0 = 



1 0 
0 0 
0 1 



-1 0 

0 p~ N 



2w = di*g{p\p\...,p»-»} t p=l- x ( 7/ ) 
. K 0 = 0, x 0 |o = x 0> Kt = p' i K k (31) 



here. 
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y k : the observation signal, 

F k : the dynamics of the system, 

H k : the observation matrix, 

x~ k | k : the estimated value of the state x k at the time k using the 
observation signals y 0 to y k , 

K s , k : the filter gain, obtained from the gain matrix K" k , and 
R e , k , lT k : an auxiliary variable, 

a step at which the processing section stores an estimated 
value of the state x k by the hyper H*, filter into the storage section; 

a step at which the processing section calculates an existence 
condition based on the upper limit value Yf and the forgetting factor 
p by the obtained observation matrix H ± or the observation matrix 
Hi and the filter gain K S/i , and 

a step at which the processing section sets the upper limit 
value to be small within a range where the existence condition is 
satisfied at each time and stores the value into the storage section, 
by decreasing the upper limit value Yf and repeating the step of 
executing the hyper H«> filter. 

12. (Canceled) 

13. (Amended) The system estimation method according to 
claim 7 or 9 or 11, wherein an estimated value z v k | k of the output 
signal is obtained from the state estimated value x~ k | k at the time 
k by a following expression: 

z\|k - H k x~ k | k . 

14. (Amended) The system estimation method according to 
claim 7 or 9 or 11, wherein the filter equation is applied to 
obtain the state estimated value x~ k | k =[h~i[k] , • • • , h~ N [k]] T 

a pseudo-echo is estimated by a following expression: 



dk = £M*Hbft k = 0,1,2, • 

t=0 



(34) 



and 

an echo canceller is realized by canceling an actual echo 
by the obtained pseudo-echo. 

15. ( Amended ) A system estimation program for causing a 

computer to make state estimation robust and to optimize a 

forgetting factor p simultaneously in an estimation algorithm, in 
which 

for a state space model expressed by following expressions: 
x k+1 - F k x k + G k w k 
y k - H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 

which indicates a ratio of a filter error to a disturbance including 

the system noise w k and the observation noise v k and is weighted 

with the forgetting factor p is suppressed to be smaller than a 

term corresponding to a previously given upper limit value Yf , and 
the system estimation program causes the computer to execute : 
a step at which a processing section inputs the upper limit 

value Yf. the observation signal y k as an input of a filter and a 
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value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value Yf; 

a step of executing a hyper H« filter at which the processing 
section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s , k by using the forgetting factor p and a gain matrix 
K~ k and by following expressions: 



K sJc = piK k (: } 1)/JUU) 







0 


0 







0 

Kk 



— 1 W ~ 

R eJc C k+lL k 



ReMl = R e>k ~ Ck+lLkRr^LkC^ 



Rr tk +l = Rr fk - L k C k + x R~ k Ck+lLk 

Where , 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



, H k+l = [u k+1 u(k + l-N)] = [u(k + l)u k ] l H, = [«(1),0,... ) 0] 



R t ,\ = Rx + C l S lVs C\ 



' R -l p 0 



L 0 = 



1 0 
0 0 
0 1 



e 7?(Af+l)x2 R _ I -1 0 



, K 0 = 0, x 0 | 0 = x 0 , Kt-p~iK t (31) 
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here, 

yit: the observation signal, 
F k : the dynamics of the system, 
H k : the observation matrix, 

x'kifc: the estimated value of the state x k at the time k using the 
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observation signals yo to y k , 

K s#k : the filter gain, obtained from the gain matrix K~ k , and 
R e ,k, L\: an auxiliary variable. 

a step at which the processing section stores an estimated 
value of the state x k by the hyper H«> filter into the storage section; 

a step at which the processing section calculates an existence 
condition based on the upper limit value y f and the forgetting factor 
p by the obtained observation matrix Hi or the observation matrix 
Hi and the filter gain K S/i# and 

a step at which the processing section sets the upper limit 
value to be small within a range where the existence condition is 
satisfied at each time and stores the value into the storage section, 
by decreasing the upper limit value Yf and repeating the step of 
executing the hyper H* filter. 

16. (Amended) — A computer readable recording medium 
recording a system estimation program for causing a computer to 
make state estimation robust and to optimize a forgetting factor 
p simultaneously in an estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+ i ■ F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k ; a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 
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which indicates a ratio of a filter error to a disturbance including 
the system noise w k and the observation noise v* and is weighted 
with the forgetting factor p is suppressed to be smaller than a 
term corresponding to a previously given upper limit value Yf / and 

the computer readable recording medium recording the system 
estimation program causes the computer to execute: 

a step at which a processing section inputs the upper limit 
value Yf* the observation signal y k as an input of a filter and a 
value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value y f ; 

a step of executing a hyper H«> filter at which the processing 
section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s , k by using the forgetting factor p and a gain matrix 
K' k and by following expressions: 
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**|* = &k-i|*-i + K s>k (y k - ffjtXik_i| fc _i) 



Kk+i 




0 


0 







0 



Re,k+1 - R e,k ~ C k +iL k R r j.L k C k+l 

- T " T _j u 
flr,fc+l = -Rr,fc ~ L k C k+l R ek Ck+\L k 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



Where, 
Cm ~ 



=ie, + C 1 I7 1 , 0 C 11 R x = 



1 0 
0 0 
0 1 



if Hi = ki ti(* + 1 - N)] = (u(/: + 1) Uk ) t H\ = [u(l),0, ... ,0] 

, i7i, 0 = diag{p 2 J p 3 1 ... 1 p^+ 2 }, p=l- x ( 7/ ) 

#0 = 0, i 0|0 = x 0 , Kk = p~*K k (31) 



P 0 
PT/J 



0 2 



-1 0 
0 p' N 



here, 

y k : the observation signal, 
F k : the dynamics of the system, 
H k : the observation matrix, 

x~ k | k : the estimated value of the state x k at the time k using the 
observation signals y 0 to y k , 

K s , k : the filter gain, obtained from the gain matrix K~ k , and 
R e , k , L\: an auxiliary variable. 

a step at which the processing section stores an estimated 
value of the state x k by the hyper H«, filter into the storage section; 

a step at which the processing section calculates an existence 
condition based on the upper limit value Yf and the forgetting factor 
p by the obtained observation matrix H± or the observation matrix 
H ± and the filter gain K S/i , and 

a step at which the processing section sets the upper limit 
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value to be small within a range where the existence condition is 
satisfied at each time and stores the value into the storage section, 
by decreasing the upper limit value Yf and repeating the step of 
executing the hyper H* filter. 

5 

I 17. (Amended) A system estimation device for making state 

/ estimation robust and optimizing a forgetting factor p 

simultaneously in an estimation algorithm, in which 

for a state space model expressed by following expressions: 
10 x k+1 = F k x k + G k w k 

y* = H k x k + v k 

z k = H k x k 
here, 

x k : a state vector or simply a state, 
15 w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 
20 G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy gain 

which indicates a ratio of a filter error to a disturbance including 

the system noise w k and the observation noise v k and is weighted 

with the forgetting factor p is suppressed to be smaller than a 
25 term corresponding to a previously given upper limit value Yf / and 
the system estimation device comprises: 

a processing section to execute the estimation algorithm; 

and 

a storage section to which reading and/or writing is performed 
30 by the processing section and which stores respective observed 
values, set values, and estimated values relevant to the state space 
mo de 1 , where in , 
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a means at which the processing section inputs the upper limit 
value y f , the observation signal y k as an input of a filter and a 
value including an observation matrix H k from the storage section 
or an input section; 

a means at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value y f ; 

a means of executing a hyper Hoo filter at which the processing 
section reads out an initial value or a value including the 
observation matrix H k at a time from the storage section and obtains 
a filter gain K s , k by using the forgetting factor p and a gain matrix 
K' k and by following expressions: 







0 


0 







0 



Rfjc+l 

Where, 



= Rr,k ~ L k C k+l R~ k Ck+lL k 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



H k +l 



. **+i as |tt*+iti(t + l-iV)] = lH(* + l)u»j l H x = [«(!), 0.....0] 



P 0 
PT/1 



0 2 



L 0 = 



1 0 
0 0 
0 1 



-1 0 
0 p~ N 



, £i|o = diag{pV,... p = l-xhf) 
, K 0 = 0, x 0 \o = x o , K t =p' i K l (31) 



here, 

15 y k : the observation signal, 

F k : the dynamics of the system, 
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H k : the observation matrix, 

x~ k | k : the estimated value of the state x k at the time k using the 
observation signals y 0 to y k# 

K s#k : the filter gain, obtained from the gain matrix K" k , and 
5 R e ,k* L~ k : an auxiliary variable. 

a means at which the processing section stores an estimated 
value of the state x k by the hyper Hoo filter into the storage section; 

a means at which the processing section calculates an 
existence condition based on the upper limit value Yf and 
10 forgetting factor p by the obtained observation matrix H± or the 
observation matrix H ± and the filter gain K Sfl# and 

a means at which the processing section sets the upper limit 
value to be small within a range where the existence condition is 
satisfied at each time and stores the value into the storage section , 
15 by decreasing the upper limit value Yf and repeating the means of 
executing the hyper Hoo filter. 
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